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We study slow-roll inflation on a three-brane in a five-dimensional bulk where the effects of energy 
loss from the brane due to graviton emission is included in a self-consistent manner. We explicitly 
derive the form of the energy loss term due to inflaton-to-graviton scattering and thus determine 
the precise dynamics of the two resulting inflationary solutions. What is also remarkable is that 
nonconservation of energy on the brane causes the curvature perturbation to not be conserved on 
superhorizon scales even for the purely adiabatic perturbations produced in single-field inflation. 
Thus the standard method of calculating the power spectrum of inflaton fluctuations at Hubble 
exit and equating it to the power spectrum at horizon reentry no longer holds. The superhorizon 
1 evolution of the perturbations must be tracked from horizon exit through to when the modes reenter 

the horizon for the late time power spectrum to be calculated. We develop the methodology to do 
this in this paper as well. 
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Recent progress in superstring theory and M- Theory has led to increased interest in cosmological models with extra 
dimensions. Braneworld models are a particular class of these models in which our four-dimensional Universe lies 
Q_j on a three-brane embedded in a higher-dimensional (possibly large) spacetime (referred to as the bulk) Matter 
fields are confined to the brane while gravity is free to propagate in the bulk. An example of these models which has 
been very well studied is the one proposed by Randall and Sundrum where a ^-symmetric (i.e. mirror symmetric) 
three-brane is embedded in five-dimensional Anti-de Sitter ( AdS) spacetime Q ■ 

However, as mentioned, gravity is not confined to the brane, and so fluctuations in the matter fields on the brane 
can emit gravitational waves (or bulk gravitons) into the bulk. The brane energy-momentum tensor is therefore not 
conserved and there will be a flow of energy from the brane to the bulk. It is essential that this effect be considered for 
a cosmological model to be realistic. And in the context of the Randall-Sundrum II model, it has been realised that 
£SJ ! the bulk cannot simply be pure AdS but must also have a bulk energy-momentum tensor describing the radiation of 
j" ■ gravitons from the brane. This has been studied by several authors. For example in [3[, the authors consider thermal 
emission of gravitons from the brane. In 0, [E] on the other hand, an exact bulk solution is used, the five-dimensional 
generalisation of Vaidya's metric, which describes gravitons being emitted orthogonally from the brane. 

These two effects, namely the presence of extra dimensions and the nonconservation of energy-momentum on the 
brane, may significantly affect our picture of the early Universe. As such, inflation (and the spectrum of perturbations 
generated) offers a powerful method of testing such a model. The spectrum of inflationary perturbations is likely to 
be modified due to these two effects. Common observables like the spectral index for this model can therefore be 
calculated and confronted with astronomical observations. 

Inflation on a brane within a higher dimensional spacetime has already been extensively studied @, [h|. However, 
to the best of our knowledge, the effect of the emission of gravitons from the brane (i.e. the nonconservation of brane 
energy- momentum) on inflation has only recently been considered [7| . The purpose of this paper is therefore to build 
on the work done in [?J and to further investigate novel consequences of nonconservation of brane energy-momentum 
on a canonical model of slow-roll chaotic inflation on the brane. It is assumed that the inflaton is confined to the brane 
and that it is the dominant source of energy-momentum. The effect on the brane-bulk system of graviton emission 
from the brane is included in a self-consistent manner, as was done in 0]. 

This paper is organised as follows. In section 2, we present the governing equations of the model. The slow-roll 
inflationary background equations are then the focus of section 3 and the precise dynamics of the two resulting 
inflationary solutions are derived. In section 4, we consider scalar perturbations about the homogeneous background. 
In particular, we demonstrate that the curvature perturbation is not conserved on superhorizon scales even for 
adiabatic perturbations and that new tools must be developed to calculate the late time power spectrum. We end 
with our conclusions in section 5. 
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II. BASIC MODEL 



A. Governing Equations 



We follow the model presented in [4(. This model consists of a ^-symmetric three-brane, with tension A, embedded 
in a five-dimensional bulk spacetime. Ordinary cosmological matter is present on the brane with energy density p 
and pressure p (later, when we consider inflation, the inflaton field will be the dominant component on the brane). 
The five-dimensional bulk spacetime has negative cosmological constant A = — 6/1 2 and an energy-momentum tensor 
given by: 

T AB = ak A k B , (1) 

where k A is a null vector orthogonal to the brane and a represents the energy flux of the bulk gravitons from the 
perspective of a brane observer. We assume that the five-dimensional Einstein equations hold in the bulk, thus the 
bulk spacetime satisfies 

G ab + A9ab = 4Tab- (2) 
Here n 2 = M^ 3 is the five-dimensional gravitational coupling. 



These five-dimensional Einstein equations can be solved using the following metric ansatz 

ds 2 = -f(r, v)dv 2 + 2drdv + r 2 dx 2 , f(r, v) = p 2 r 2 - (3) 

This is the five-dimensional generalisation of Vaidya's metric [l2| which describes the spherically-symmetric spacetime 
surrounding a radiating star. The coordinate r is identified with the scale factor a on the brane. ft should be noted 
that the price we have to pay for having an exact analytic solution for the bulk spacetime is that we must assume that 
all the gravitons are emitted "radially" from the brane, an assumption which is not likely to be true at high energies. 
(Here, the radial direction r essentially means the direction orthogonal to the brane in the bulk spacetime.) 

Using various components of the five dimensional Einstein equations, as well as the Israel junction conditions at 
the brane location, we obtain the following equations: 



H 2 = 4 2 + 4 x+ c^ (4) 

36^ 18 H a 4 v ' 



- -§-a 3 {a- Vf + a 2 ) (5) 



3H(p + p) = -2(70 (6) 



Before we proceed, some comments should be made regarding this system of equations. First, for the modified 
Friedmann equation Q, if we use the Randall-Sundrum condition (k 2 X = 6p) and make the following identification, 

^ = n 2 p = SttG = M~\ (7) 

then the second term agrees with the standard Friedmann equation. Also, by comparing the magnitudes of the first 
and second terms, we can define two regimes: the high energy regime (p 3> A) where H 2 oc p 2 and the low energy 
regime (p A) where H 2 cx p. The radiation-like term originates from the projection of the bulk Weyl tensor 
and hence C(v) is often referred to as the Weyl parameter. If the bulk is empty, then C(v) is constant and the 
5-dimensional metric is equivalent to the AdS-Schwarzschild metric (C will therefore be a measure of the mass of the 
bulk black hole). Here, we have also chosen to set the four-dimensional cosmological constant to zero; we do not 
attempt to solve the cosmological constant problem with this model. 



The vector k A and the unit normal n A which would appear in the junction conditions have been selected and normalised appropriately 
so that these equations appear in the above nice form. More details on this are given in Q. 
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And lastly, for the nonconservation of energy equation ([6]) , the term on the RHS represents the loss of energy from 
the brane into the bulk. The factor 2 is present because the brane is radiating on both sides, into two copies of the 
same bulk spacetime. 

We are primarily interested in a period of slow-roll inflation during the high energy regime. In this limit, the 
previous system of equations can be solved. We find that the term is just a very small modification to the 
Friedmann equation during slow-roll inflation. We therefore choose to neglect it. The system of equations reduces to 
just 

H 2 = |p 2 (8) 
p + 3H(p + p) = -2d. (9) 



B. Calculation of Loss Term due to Graviton Emission 



We now need to determine the energy loss a. We assume that all the energy loss from the brane is due to the 
process <fi + 4> —¥ G, i.e. the emission of a bulk graviton from the scattering between two inflaton particles. Our 
calculation closely follows section III in reference [4j. We make several simplifying assumptions. First, it is assumed 
that the cosmological expansion can be neglected in the computation of the cross-section for this process. A Minkowski 
background in the brane is therefore assumed. The bulk gravitons can thus be described as (generalised) Kaluza-Klein 
modes (as done in Q). The required amplitude for the process, <j) + (f> — > G, is then given by [1, H| 



\M\ 2 = — K ls 2 (10) 
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where the Mandelstam variable is s = — (E^ + Ek 2 ) 2 + (ki + k2) 2 . ki^ are the momenta of the two incoming inflaton 
particles and E kl 2 represent the energies of the inflaton particles. 

To determine the energy loss due to this process, we need to use the relativistic Boltzmann equation with a collision 
term. However, since we are interested in the energy density on the brane and not the number density of inflaton 
particles, we must multiply the collision term by the graviton energy and integrate over all Kaluza-Klein modes and 
over momentum space. The following energy nonconservation equation is obtained: 

d 3 k r 

P 



where 



3H(P + P) = -Ijd m d j^c (11) 

C = J {2n)3 2 Eki {2 ^2E k2 l^lV^( 2 ^)^ (3) (k m k i k 2 )(2^(ff km - E ki - E ka ). (12) 

Here, k m is the bulk graviton momentum and n kl 2 are the number densities of inflaton particles with momenta kx 2 - 
Similarly, the graviton energy is given by £k m = -v/lkml 2 + m 2 where m which labels the Kaluza-Klein mode can be 
thought of as the graviton mass from the brane perspective. The inflaton particles are treated as effectively massless, 
so E kl 2 = |k 1)2 |. 

During inflation we assume that the inflaton fluctuations are in the standard de Sitter vacuum, i.e. the Bunch- 
Davies vacuum. We further note that for a scalar field in de Sitter space, any timelike observer in the Bunch-Davies 
vacuum state will measure a thermal distribution of particles with temperature, T = ^ @. We therefore take the 
number density of inflaton particles to be the thermal distribution: 

= -ST- • (13) 



The RHS of (jTTj) can then be integrated explicitlj0 to give 



7 . 



2 To obtain this explicit result we must further assume that the contribution to the energy loss is dominated by the "heavy" gravitons, 
i.e. those with m 2> /«■ This is explained in reference Q. 
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a, defined implicitly above, is approximately 3 x 10 . The p° term is also likely to be very small in magnitude. 
However, k^ 8 which is a measure of the difference between the 4D Planck mass, M p and the 5D fundamental mass 

scale, M 5 = /% 2 ^ 3 can potentially be quite large. For example, taking V — ^m 2 (j) 2 and M 5 = 10 _3 M p , we find that 

a4 8 p 8 .( Mr ^ 8 



Thus, the loss term will generally be small, but not necessarily to such a negligible extreme. The dependence of K5 8 
on the relative sizes of M5 and M p has been included explicitly in equation (fT5j). 



III. SLOW-ROLL INFLATION ON THE BRANE 
A. Two Different Branches of Slow-roll Dynamics 

We consider an inflaton field confined to the three-brane. It is assumed that the inflaton field <j> (with self-interaction 
potential V((f>)) forms the dominant contribution to the brane energy density and pressure. Our work in this section 
builds on ■ In 0] , the authors discovered that energy loss will result in the presence of a new inflationary branch. 
Here, we explicitly derive the dynamics of both branches, examine the stability of the new branch, and investigate 
why a new branch appears when there is energy loss. 

We first observe that since the nonconservation of energy equation and the modified Friedmann equation are different 
to those of standard 4D general relativity, the condition for inflation, i.e. the condition to have accelerating expansion, 
is now different as well. The following inequality is obtained 

a > 0^>p < ~\p(l + ianfp G ) , (16) 

which is a stronger condition than the standard 4D result p < — |p. If we assume there is no energy loss (a = 0), 
then we obtain p < — |p which is in agreement with the result obtained in Q for the high energy regime. 
With the usual expressions for the energy density and pressure of the inflaton scalar field, 

P=\tf + V (17) 

P = \<P 2 - V, (18) 

the nonconservation of energy equation becomes: 

4>$ + 3Hj> + V) = -aKf{^ 2 + Vf. (19) 

To obtain a negative pressure fluid, we employ the slow-roll mechanism. It is assumed that the <j> and ^<j) 2 terms are 
negligible compared to the remaining terms. We thus obtain the following equation for <f> 

3H4> 2 + V'4> + a K l 8 V 8 = 0. (20) 

This quadratic equation in can be solved to give 



(f, = -6H ± mi 1 P— (21) 

Due to the smallness of the loss term, the second term in the squareroot is much smaller than the first. The squareroot 
is therefore expanded to O(a) to give the two following solutions, 

l--V\°tfV (22) 



3H V 

which essentially is just the (slightly modified) standard slow-roll dynamics, and 



(23) 
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which we dub the nonstandard inflationary branch. 

So it is clear that the emission of gravitons (which causes the brane energy-momentum to not be conserved) has 
led to the possibility of having two different branches of inflationary dynamics: the modified slow-roll branch and the 
nonstandard inflationary branch. 



B. Is the Nonstandard Inflationary Branch a Stable Attractor Solution? 

We would like to determine whether the nonstandard branch is stable to small perturbations. If it is, then we 
expect there to be a basin of attraction for this branch. The dynamics of the nonstandard inflationary branch should 
therefore be observed for a range of initial conditions. Or, in other words, the nonstandard branch will thus be a 
stable attractor solution. 

To do this, we consider an inflationary solution starting slightly away from the nonstandard branch, <j> + and check 
analytically whether the perturbed solution converges to the nonstandard branch. For the purpose of this calculation, 
we specify to potentials of the form, 



Writing <fr as 
the equation for 

becomes 



V = V(4>) = H n - (24) 
= 0++ <50, (25) 

■3H<j> + V') = -anl 8 (^<j) 2 + V) S (26) 

(0+ + 50)(0+ + 60 + 3tf (0+ + 8cj>) + V) = -anfC-ih + H? + Vf (27) 

We neglect second order terms, <50 2 ,84>8(f>, etc. Also, we are only interested in the solution up till the end of inflation, 
so we assume that the slow-roll conditions hold for the nonstandard branch. We further define V+ = V(<f>+) = /30™ 

and H + = ^ (|0+ + V+j . The previous equation then reduces to: 

+ c50 + 3#0 + + ?>H4> + 8<j> + <p+V + 3HS<t><f> + + 8<f>V = -anf{^ + 8<t)V 7 + V 8 ) (28) 
We now need to determine V,V, and H, 

V = /?0" = /3(0+ + ,50)" = /3<f> n + + Prupf^St = V+ + pri^ 1 8<j> (29) 

V = V[ + n(n - l)/30!p 2 (50 (30) 
H = f + V) = f ( 2^ + + ^ + V) = f ( 2^+ + ^ + V++ 

2 2 

= H + + ^np^H + y0+<50- (31) 
Substituting these into the previous equation and neglecting second order terms gives: 

2 2 

4>+5<f> + 3H + ^l + ^nP^ft+H + y 0+<50 + 6H+4>+S<f> + 0+ V + + 0+n(n - l)/30Ip 2 c50 + <W+ 

= -anf (84>+6<pVl + Vf + 8V_£V|<ty) . (32) 

For the nonstandard branch, we have that 3£f+0+ + V+0+ + an\ & Vf_ = and that + = Using the first of 

these equations, the equation for 8(f) reduces further to: 

K 2 

0+60 + -^■n/30Ip 1 0+60 + k\V + 4> + 84> + 4>+n{n - l)Pcf>\- 2 8cj> + 5(f>V[ 
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= -anf (8<p+S(j)Vl + 8V^Vj(5</>) . (33) 
Using V + = P4>+ and V\_ = nj3(j> , ] r ~ 1 , and after some algebra, we obtain: 



+ + n(n IWV 2 + 8a ^ n ^ Sr j 54> = 0. (34) 

Substituting for <p + finally gives us: 

H + i^lm + aKl*0 7 <fr - — 15^+2 ) ^ + (- ^Y^ S <P+ -n(7n+ 1)^) H = 0. (35) 

We now need to solve this 2nd order ODE to determine the behaviour of 5(f). Since we are only interested in whether 
<f> converges to (or diverges from) <j> + , we assume the slow-roll conditions for <p + hold (otherwise <fi + would no longer 
be a valid solution). With this in mind, we treat <f) + as roughly constant, and so we just have an ODE with constant 
coefficients. 

And so we have solutions of the form e At where A is either negative or very close to zero. We therefore see that 
we have a decaying solution and another that is almost constant (which just corresponds to a shift along the inflaton 
trajectory). Hence, the nonstandard inflationary branch is a stable attractor solution. 

C. Why is there a Second Branch? 

In this subsection we consider the origin of the second inflationary solution. It is straightforward to understand why 
the modified slow-roll branch exists; the energy loss term perturbs the equation of motion for <j> slightly which then 
slightly perturbs this inflationary solution away from the standard slow-roll dynamics. But why is a second solution 
introduced? 

In actual fact, the canonical equation of motion for the inflaton (i.e. with no energy loss) 

</> (J> + 3ii"0 + V") = 0, (36) 

has two solutions, 

4> + 3H<j) + V = 0, (37) 

which reduces to the standard slow-roll dynamics if we assume the slow-roll conditions hold, and 

4> = 0. (38) 

This second solution is typically neglected because if 4> is constant then there is no dynamical way of ending inflation 
and the universe would stay in a de Sitter phase indefinitely. However, the presence of energy loss perturbs this 
solution away from a constant <j> solution, to give: 

„...18t/8 

(39) 



V 

Thus, with energy loss, a second inflationary solution is allowed even for the canonical inflaton action. 

IV. CALCULATION OF THE SPECTRUM OF SCALAR PERTURBATIONS 
A. Curvature Perturbation Not Conserved 



The standard way of relating the fluctuations during inflation to late time observables is to calculate the power 
spectrum of quantum fluctuations at horizon exit and to conclude that this is equal to the power spectrum at horizon 
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reentry at late times if we assume that the curvature perturbation is conserved on superhorizon scales (for the adiabatic 
perturbations which are generated in single-field slow roll inflation). 

However, this is no longer the case when energy is not conserved on the brane (due to graviton emission, for 
example) and we therefore need to resort to different methods. The evolution of the perturbations must be tracked 
on superhorizon scales until horizon reentry at late times. In this subsection, we demonstrate that the curvature 
perturbation is not conserved on superhorizon scales even for the adiabatic perturbations generated in single-field 
slow-roll inflation. In the following subsections, we go through in detail the method we have used to calculate the late 
time power spectrum. 

To demonstrate that the curvature perturbation is not conserved, we work in the Newtonian gauge for simplicity. 
Thus we use the following form of the perturbed Friedmann metric on the brane: 

ds 2 = -(1 + 2^)dt 2 + a 2 [(l - 2#)£y]dxW. (40) 

Only the scalar perturbations are of interest to us for this analysis, so we neglect vector and tensor modes. For the 
calculation of the time derivative of the curvature perturbation, we require the nonconservation of energy equation. 
Starting from the nonconservation of energy-momentum on the brane, we act with u" , a timelike vector of a comoving 
observer on the brane, to obtain nonconservation of energy: 

u v V^ v = -anfp 8 . (41) 

And if we are interested in just the first order perturbation terms on superhorizon scales, we then hav^l 

8p + 3H(5p + Sp) - 3(p + p)# = -8anl 8 Spp 7 - Qanfp 8 . (42) 

The last term on the right arises from the perturbation to u in the term u v V \ l T^ 1 ' '. 

The zeroth order equation (which gives the background nonconservation equation) will also be useful: 

p + 3H(p + p) = -cweJV- (43) 

We want to calculate the time derivative of £ which is the gauge-invariant variable representing the curvature pertur- 
bation on uniform density hypersurfaces, £ = — \f r — ^-Sp. 



HSp anfp 8 8p 



m{p + p) + anfp% 6(3H(p + p)+aK 1 5 8 p s ) 
3H 2 (l + |) Sp ZaHnffSp 



3H(p + p) + aK 1 5 8 p 8 3H(p + p) + aK 1 5 8 p 8 ' 
The nonconservation of energy equation can also be rewritten as: 

H5p 3H 2 Sp 3H 2 Sp 



(44) 



m{p + p) + anfp 8 3H(p + p)+anl 8 p 8 ZH (p + p) + anf p 8 
3H(p + p)^ %an\ 8 H8pp 7 ^Hanfp 8 



3H(p + p) +aKl 8 p 8 3H{p + p)+aK 1 5 8 p s 3H{p + p) + anfp* ' 
We require the definition of the non-adiabatic pressure perturbation, Sp na d, as well, 



(45) 



Spnad = Sp-^-. (46) 
P 



We can now combine the previous three equations to obtain: 

i = 3H 2 5p nad anfp 8 ^ a^K 20 p 9 

4 3H(p + p) + anfp 8 2,H{p + p)+an\ 8 p 8 3H (p + p) + anf 'p 8 ' 
And for adiabatic perturbations on superhorizon scales, the first term vanishes, leaving: 

a.K\ 8 p 8 i! a^nfp 9 
3H(p + p)+aKl 8 p s ~ 3H(p + p)+a K l 8 p 

Thus, we see that the curvature perturbation is not conserved on large scales, even for the purely adiabatic pertur- 
bations generated in single field inflation. 



We are now working in momentum space, but for notational convenience, the fc-dependence has been suppressed. 
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B. Quantisation of Inflaton Fluctuations 

In the previous subsection, we have demonstrated that the curvature perturbation is not conserved on superhorizon 
scales. The power spectrum at horizon reentry at late times therefore cannot merely just be matched to the power 
spectrum of quantum fluctuations at Hubble exit. New methods are required to track the superhorizon evolution of the 
curvature perturbation through till horizon reentry during the late time low energy regime. To achieve this, we adopt 
a method commonly used in multifield slow-roll inflation models (where here the non-adiabatic perturbations cause 
the curvature perturbation to evolve on superhorizon scales). Essentially, the method is as follows: the perturbation 
equations are solved on superhorizon scales, then, to fix the integration constants that arise, the amplitude of the 
quantum fluctuations at horizon exit is used as an initial condition. We now have a solution to the perturbation 
equations that is valid on superHubble scales, so we can use this to evaluate the power spectrum at horizon reentry 
at late times. 

In actual fact, our case is slightly more complicated than this due to the presence of the two different regimes: the 
high energy regime and the low energy regime. More specifically, the perturbation equations must first be solved on 
superhorizon scales in the high energy regime under slow-roll inflationary conditions and assuming that the adiabatic 
mode is weakly time-varying. The integration constants which arise arc then fixed by comparison to the power 
spectrum of quantum fluctuations at Hubble exit. We assume that this period of high energy slow- roll inflation 
is followed by a low energy radiation era in which standard 4D general relativity holds and standard cosmology is 
recovered. And once we have recovered standard 4D general relativity, we know that the curvature perturbation 
£ will be conserved on superhorizon scales (for the adiabatic perturbations in single- field inflation) [l3l . [l8| . So the 
superhorizon solution found previously which is valid only during the high energy slow-roll inflationary period can 
be used to calculate the value of £ at the end of inflation, which will then be conserved until horizon reentry at late 
times. The power spectrum of £ at late time horizon reentry can thus easily be obtained. 

In principle, when performing this calculation, the full system of bulk and brane metric perturbations together 
with the inflaton perturbation should be considered in its entirety. We cannot just perturb about a 4D Friedmann 
metric but should include the coupling to the bulk metric perturbation as well. This complicates our calculation 
in two ways. First, to obtain the power spectrum of scalar perturbations at horizon exit, we should quantise the 
complete set of perturbations for coupled brane-bulk system and not merely perturb about a 4D Friedmann metric 
and neglect bulk effects. However, it was shown in 10] that even when the coupled system of brane-bulk scalar 
perturbations is consistently quantised, the bulk metric perturbation only introduces corrections at first order in the 
slow-roll parameters. So, in this instance, we will neglect the bulk metric perturbation, and quantise the system of 
inflaton-brane metric perturbations only. Admittedly we are using a different bulk metric to [10( , but we assume that 
any additional corrections introduced by our bulk metric will be further suppressed by a, the coefficient of our loss 
term, and this simplifying assumption should therefore remain valid. 

The second way in which the coupling to the bulk metric perturbation affects our calculation is its effect on the 
superhorizon evolution of the brane scalar perturbations during the high energy inflationary regime. When considering 
how the scalar perturbations on the brane evolve, we should again consider the full brane-bulk system of equations. 
However, it was shown in [l3j that for any braneworld model where the 4D energy-momentum tensor is conserved, 
then the brane curvature perturbation is conserved on superhorizon scales (for pure adiabatic brane perturbations). 
So for our model, we assume that it is the nonconservation of energy on the brane and not the higher-dimensional 
bulk perturbations which is the dominant effect on the superhorizon evolution of the brane scalar perturbations. We 
therefore neglect the effect of the bulk perturbations for the rest of this section. 

In this subsection, we aim to quantise the brane scalar fluctuations to determine the power spectrum at Hubble exit. 
There are two ways in which we could proceed. First, we could start from the second order action, change variable 
appropriately, and quantise. However, in our case, while we could start from a 5D action, we have made several 
simplifying assumptions regarding the bulk effects, so perhaps this is not the best way to proceed. Alternatively, we 
could start from the scalar perturbation equations, define a Mukhanov variable [l5j such that we obtain the harmonic 
oscillator equation, and then quantise. By right, since the equations we obtain are linear, any new variable defined 
through some numerical rescaling would be a viable solution as well. Thus, to ensure that we have chosen the right 
normalisation, we demand that, in the a — > limit (where standard inflationary theory is recovered), the standard 
solution to the mode equation is reobtained. 

We first derive the perturbation equation for <5</@. We start with the nonconservation of energy equation, 

u v V^ = -an^p 8 . (49) 



4 Our method in this subsection is similar to 
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and then focus on the first order perturbation terms which gives the equation 



k 



Sp + 3H(5p + Sp) - 3$(p + p) j Sq = -8anl 8 6pp 7 - $an 



18-8 
5 P 5 



(50) 



in the Newtonian gauge and for zero anisotropic stress (which is true if we neglect bulk effects). 
Using the expressions for Sp, Sp, and Sq in slow-roll inflation 

Sp = 4>5cf> - $</> 2 + V'S(j) 



Sp = 



$0 2 - V'S(f>, 



(51) 
(52) 



Sq = - 



(53) 



the perturbation equation becomes 



4$0 2 - 2$<j)4> + V"4>5(t> + V'5<t> + 6H<])S<j) - 6£T$0 2 + 



= — 8a He. 



<i>/;- - \ "s, » ( x -i\? + v) - $«4 8 + v" 



(54) 



This equation can be simplified slightly by using the background equation, 



+ 3H(j) + V') = -anf ( -</> 2 + V 



(55) 



to remove the second time derivatives of <fi. It simplifies to 



SHStf) - + V'Scj) + 2V'<f> 



4 8 ( \i> 2 + v 



13 • 3$U 15-, 84>V 8V'i 

$0 : S(b + -X - 

2 Y 6 2 Y d? 6 



(56) 



which is our starting point. 

In the extreme slow-roll limit, all mass- and ^-dependent terms (which are 0(e, 77)) may be neglected^. The previous 
equation reduces to: 



3HS<f> 



a4 8 V 7 



VS<f> 8V'l 



(57) 



Changing the independent variable to conformal time gives: 



2a' anfV 8 a 



= 0. 



(58) 



We now wish to solve this equation. The standard strategy is to change variable (to the Mukhanov variable) for which 
the equation changes to the canonical form of the harmonic oscillator equation, which can then easily be quantised. 
If we are given an equation of the form 



5<t>" + AScj)' 



= 0, 



(59) 



This is equivalent to using the pure de Sitter result at horizon exit when calculating the power spectrum of quantum fluctuations. 
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we introduce the Mukhanov variable, v = z56, such that the equation for v will have no first derivative (so essentially, 
we will then have the harmonic oscillator equation). The condition for this is simply: 

97' 

— = A. (60) 

z 

With this in mind, we define 

_ f „.lBy8 a 

z = ae J *P ; where v = z56. (61) 



We now have the following equation 



v" + 



k 2 + 



v = 0. (62) 



The solution to this should be proportional to the solution of the mode equation for canonical inflation, up to 0(a) 
corrections, that is: 



e ik n , l 



" x 7sK 1 si +<Ka) - (63) 



The normalisation of this solution is therefore still arbitrary. We fix the normalisation by demanding that the 
canonically normalised solution is obtained in the a — > limit. Thus, 



e ikr) , l 



_ . 1 - — + 0(a). (64) 

Now it is not easy to write down an exact analytic expression for the O(a) term, so to make further progress 
analytically, we use the approximate a = solution for v (acknowledging that this will introduce an error into our 
final result). 

Following the standard procedure, we obtain the following power spectrum (evaluated at horizon exit) 

(66 k 66 k ,) = (2n)H(k + k')^ 

= (2n) 2 S(k + k')^eJ (65) 
So, the dimensionless and dimcnsionful power spectra (A^and Pg^ respectively) are given by 

Surprisingly, there is an explicit dependence on scale, k. 



C. Solution of $ During the High Energy Inflationary Period 



We now wish to solve the perturbation equations during the high energy inflationary regime on supcrhorizon scales. 
The perturbation equation for 6<j> on superhorizon scales is given by, 



3HScj) - 4$0 + V'' 



k 2 



an 



18 



13 • 3$V 15-, 5$V 8V54) 

$0 : 56 + —r r-t- 



(67) 



We want to solve this equation under slow-roll conditions and assuming that the adiabatic mode is weakly time- 
dependent so that we can neglect the time derivative of $ as well as second time derivatives. The equation simplifies 
further to 



3HS6 + V"5cf> + 2V'$ = anfV 1 



56V Z<S>V Wt 



(68) 
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Note that unlike in the previous subsection, not all the ^-dependent terms disappear, hence we combine this equation 
with the Oi-componcnt of the Einstein equations to obtain a closed set of equations. 
The Oi-component of the Einstein equations is given by d 



+ $ = 

which under the assumption that $ is negligible reduces to 



(69) 



(70) 



We solve this coupled set of equations using the method in [Iff . We first substitute the Einstein Oi-equation into the 
first order perturbation equation of energy nonconservation, which results in: 



V" 



2V'd> 3an 18 V 8 8an 18 V 7 V' 



H H 

The independent variable is then changed from time to <fi (that is, 



5<f> = 0. 



(71) 



= 8<b' 6 where / denotes differentiation wrt to 



8(f>- 



V" + 



2V'd> 3an\ 8 V 8 8an 18 V 7 V 



H 



H 



0. 



(72) 



Before continuing to solve these equations, we first check to see if this result agrees with standard low energy inflation 

r 

V" 2V' 



(i.e. the case with a = and H 2 = %V). Under these conditions, our previous equation reduces to just 



_3Hj> 3H 2 

And using standard inflationary dynamics 3H(f> = — V', this becomes 

v v 



= 0. 



0. 



This first order differential equation can easily be integrated to obtain 

A V' 

-, A some constant 



V 



(73) 



(74) 



(75) 



which is the correct result for standard low energy inflation [171 ] . 

We now solve the perturbation equations for our inflationary dynamics, i.e. both for the modified slow-roll branch 
and the nonstandard inflationary branch. 

The background dynamics for the modified slow-roll branch are given by 



-V anfV 8 



3H 



V 



(76) 



With this and the high energy result for the Friedmann equation (H — -|— ), the perturbation equation becomes (to 
leading order in a) 



-V 



6anl 8 V 8 H 



V 



+ 



V" 



2V 5an\ 8 V |s - 

1 5 OAn,k-^ rl 

3H 2 H 



2Aan^V'H 



0. 



Dividing through by the coefficient of the first term, we obtain 

-V" 2V 



V 



3H 2 



aL 



60 = 0, 



(77) 



(78) 



By right, there should be additional terms related to the bulk on the RHS but as mentioned earlier, we choose to neglect bulk effects 
since we assume that nonconservation of energy will be the dominant effect. 
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where aL 



cxk1 8 V 8 , 2AaK,\ s V 7 H 6aK,l 8 V s HV" 



HV T V 

Again, this first order differential equation can be integrated to obtain: 

8<j> = CV'e^e'J aL ,some constant C. 
And for the nonstandard branch, the background dynamics are given by 



V 



To leading order in a, the perturbation equation becomes 

"V" 8 V 



6' + 



V 



v 



+ aM 



0. 



This can be solved to obtain 



D— e J"*^ some const am D 



(79) 



(80) 



(81) 



(82) 



where aM — 



+ 



+ 



The constants C and D from the solutions of 8<f> on superhorizon scales in the high energy inflationary regime can 
now be determined by comparison to the power spectrum of 8<j) at horizon exit. We have the following expressions 
for 5<j) for the modified slow-roll branch and the nonstandard branch respectively 



= CV'e< v e 



-J aL 



D^e~I aM , 



Rearranging gives 



C 
D 



Me ^ e J aL 



k=aH 



8<t>V' f aM 



k=aH 



where, at horizon crossing, the quantum fluctuations of 6<j> are given by 

8<f>< 



h r a * C 



(27T) 



(83) 
(84) 

(85) 



So the high energy inflationary solutions for 8<f> on superhorizon scales for both the modified slow roll branch and the 
nonstandard branch are given respectively by: 



8cj> 
8(f> 



V'e K i v e 



— r 

V -J aL 



k=aH 



S<pV' f aM] V 8 - f aM 



k=aH 



And using the relation for $ 



H$ = 



we find the following expressions for $ in the modified slow roll and nonstandard branches respectively, 



$ = 



4>V f e K s e 



k=aH 



H 



5<bV' f aM~\ d>V 8 - f aM 



k=aH 



(86) 
(87) 

(88) 

(89) 
(90) 
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D. Matching the Solution to the Low Energy Era 



It must be stressed that these solutions are only valid on superhorizon scales during the period of high energy slow- 
roll inflation. However, we are interested in tracking these solutions until the modes reenter the Hubble radius at late 
times. As mentioned earlier, we assume that the period of high energy inflation is followed by a low energy radiation 
era in which standard 4D general relativity holds. And once standard general relativity is recovered we know that £, 
the curvature perturbation on uniform density hypersurfaces, will be conserved on superhorizon scales [l3l [l8( | . We 
therefore evaluate £ at the end of inflation (which takes place at the High Energy-Low Energy Transition (HE-LE)) 
using the superhorizon solutions found in the previous subsection. £ will then remain constant on superhorizon scales, 
therefore the expression for £ at the HE-LE transition will be valid until horizon reentry at late times. 

We will now evaluate £ at the HE-LE transition using the superhorizon solutions that we found previously. 

£ is given by the following expression 



£ = $ 



HSp 



(91) 



which during inflation can be expressed as 



$ + ^- + 0(a 2 ) 



(92) 



Using the solutions for <£> and 8<p, we obtain the following results: 



24 r. 



k=aH 



y! k 4 v — I aL 
'e s e J 







H~ 




HE-LE 


H 


0_ 


HE-LE 



— e s e 



JaL 



k=aH 



\k=aH 





H 


HE-LE 





HE-LE 



if inflation proceeded along the modified slow-roll branch, or 



£ 



j ctM~| [vf,g— / aM 



k=aH 



HE-LE 



H 



H 



HE-LE 



V' f aM~\ \ V s — f aM 

J k—aH L 

if inflation followed the nonstandard inflationary dynamics. 



HE-LE 



HE-LE 



Mk=aH 



\k=aH 



\k=aH 



(93) 



(94) 



E. Calculation of Power Spectrum 



To obtain the power spectrum of £ at horizon reentry, we simply have to calculate the power spectrum at the HE-LE 
transition since £ will be conserved until horizon reentry. Using the results of the previous subsection, we therefore 
have: 



A C = 



77 e s e 



k=aH 
1 



k=aH 



24 p 

24 p 

e & e J 



1 2 






2 






HE-LE 


'L J 
H ^ 


H 
" ? _ 


HE-LE 


A 2 


k=aH 



1 2 


H 


2 


A 2 , 


HE-LE 


A. 


HE-LE 





iH 



if inflation proceeded along the modified slow-roll branch, or 



A? 



V' f aM 



V' f aM 



k=aH 
2 



V± - J aM 



HE-LE 



£ j- K 

H d> 



k=aH 



V s - f aM 1 2 



HE-LE 



H 



2 

HE-LE 
2 r 

A 2 

HE-LE L *' 



A 2 



k=aH 



k=aH 



(95) 



(96) 



if inflation followed the nonstandard inflationary dynamics. Here, A 2 ^ is the power spectrum of inflaton fluctuations 
at Hubble exit, which we obtained earlier as 
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F. Calculation of Spectral Index 

In this section, we aim to calculate the spectral index given by the formula, 

n. - 1 



dlogA^ 



G?log k 

We first need to calculate the slow-roll parameter, e, given by 

H dlogH 
£ "~JP ~ dN ' 

Using the nonconservation of energy equation 

<j>(4> + 3H<j) + V') = -anfV*, 
and the time derivative of the Friedmann equation 



H 



we obtain the following expression, 



4£_ q K fV 8 
2H + 6H 2 ' 



(98) 



(99) 



(100) 



(101) 



(102) 



We can now use this to obtain expressions for V' for the modified slow-roll branch and the nonstandard branch, given 
by: 



and 



kI 5 



, 2 _ n 6 5 a 2 6 15 3H 15 



(103) 



(104) 



respectively. 

And substituting these into the previous expressions for the power spectra give 



18-H" 3 



5- -ra6 s 3K 2 5 H 9 



e s e 



H 2faL 



k=aH 



I- 



iHE—LE ■ 



(105) 



k=aH 



for the modified slow-roll branch and 



A/ 



^6 n 2ol6r>15 tt15 ,,32 r 
— ft J 



k=aH 



Jj2 2H^ 



(27T) 



,,.20 t -8 



iHE—LE ■ 



(106) 



k=aH 



for the nonstandard branch. The terms evaluated at the HE-LE transition have not been specified because they will 
not contribute to the spectral index as they do not have any scale-dependence. 
If we act withQ 



d 



dN d 



d log k d log k d log k dN 



(107) 



7 We use this expression since there is explicit fc-dependence as well fc-dependence implicit from the fact that different modes exit the 
horizon at different times. 
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on these expressions for the power spectra, we then obtain 



-3e + eri - ea3 9 2 7 K iH e 8e 2y/2aL , f . 74rT ^h 
"'- 1 = - 2e + e + a3 WH« 6 ^ + W (£ " a6 ^ ) 



+ (e - 6 7 cmf if 6 ) + J (e - 6 7 cmf if 6 ) ' ( j 



for the modified slow-roll branch, and 



2\/2aM , 7 4 6 . i 

n s — 1 = — e — »j H y (e - a6 k^H ) 

H 2 K5 

3 8 2 7 aniH 6 f 3 s 2 7 anikH Q 



(e-6 7 anjH e ) J (e-6 7 a K 4 5 H e ) 



(109) 



for the nonstandard branch, where i] = dloge 



dN ■ 



G. The a -»■ limit 



We find that the power spectra and spectral indices that we obtain when there is energy loss from the brane are 
quite different from the standard expressions. For standard inflation in the high energy regime, we should obtain 



for the power spectrum, and 



H 2 H 2 



k=aH 



-3e-r). 



(110) 



(111) 



for the spectral index. 

However, the power spectrum and spectral index obtained for the modified slow-roll branch should reduce to the 
standard expressions in the limit a — > 0. It is easy to see why this is true for the power spectrum. In the limit a — > 0, 
should reduce to its value at horizon exit, since if there is no energy loss from the brane, £ must be conserved on 
superhorizon scales. All terms evaluated at the HE-LE transition should therefore be able to be expressed as their 
corresponding values at horizon exit (k = aH) plus 0(a) corrections. The power spectrum for the modified slow-roll 
branch can thus be expressed as 



which is essentially just 



V 



!4 r 

W \ otL 

5 P,J 



2 

k=aH 



I 



aL 



1 2 


~H~ 


- k=aH 


J. 



k=aH 



(112) 



And this clearly reduces to 



~H~ 


2 


J. 


k=aH 







H 2 

(2n) 2 



J- 



H 2 H 2 
d> 2 (2tt) 2 



+ 0(a). 



k=aH 



(113) 



(114) 



k=aH 



if we take a — > 0. 

Note that this argument does not apply to the nonstandard branch because <j> oc a and the expression for the power 
spectrum will diverge in the limit a — > 0. 

For the spectral index however, if one works directly from the expression for the spectral index of the modified 
slow-roll branch, it is less straightforward to recover the standard expression for high energy inflation. This is because 
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energy loss, which results in superhorizon evolution beyond horizon exit, in a sense causes a loss of scale dependence 
because each mode must be evolved to and then evaluated at the same point in time (the HE-LE transition). In 
contrast, when there is no energy loss, each mode is evaluated when it crosses the horizon, and this introduces some 
terms in the expression for the spectral index which are not present for the case with energy loss. 

But if we work instead from the power spectrum, we clearly see that the standard expression is obtained, 



H 2 H 2 
02 (2tt) 2 



k=aH 



And using the expression for e for standard high energy inflation, we can reexpress this as, 

' H 3 ' 



A? 



k=aH 



We thus see that the standard expression for the spectral index is recovered, 

n s — 1 = — 3e — rj. 



(115) 



(116) 



(117) 



V. CONCLUSIONS 



We have studied slow-roll inflation on a brane where there is energy loss into the bulk due to graviton emission. 
We have built on previous work [7j and have discovered another novel result. 

First, as obtained in Q, there are two different branches of inflationary dynamics; one a slightly modified form 
of the standard slow-roll dynamics, while the other very nonstandard. In this paper, we have explicitly derived the 
form of the energy loss term from inflaton-to-graviton scattering and have obtained the precise dynamics of the two 
branches. 

Further, we find that the curvature perturbation on superhorizon scales is not conserved even for the purely 
adiabatic perturbations generated in single-field inflation. While it is common for the curvature perturbation to not 
be conserved in multifield inflationary models due to the presence of isocurvature modes, we believe this is the first 
time such a phenomenon has been demonstrated for a model of single-field inflation. Indeed, nonconservation of 
curvature perturbation on superhorizon scales should be a general feature of any single-field inflation model whenever 
the (brane) energy-momentum tensor is not conserved. 

The standard method of matching the late time power spectrum to the spectrum of quantum fluctuations at Hubble 
exit therefore no longer holds. The evolution of perturbations on superhorizon scales must be tracked until horizon 
reentry. To do this, we have adapted methods used in multifield inflationary models to our model of single-field 
inflation in a braneworld scenario. The power spectra and spectral indices obtained for the modified slow-roll branch 
and the nonstandard branch are quite different from the standard expressions for high energy inflation. But we note 
that the standard results for high energy inflation in braneworlds is recovered for the modified slow-roll branch in the 
limit a — > 0. 
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